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Abstract. An effect is considered of alternating (high-frequency) current on the 
spin-valve type magnetic junction configuration. The stability with respect to small 
fluctuations is investigated in the macrospin approximation. When the current 
frequency is close to the eigenfrequency (precession frequency) of the free layer, 
parametric resonance occurs. Both coUinear configurations, antiparallel and parallel 
ones, can become unstable under resonance conditions. The antiparallel configuration 
can become unstable under non-resonant conditions, also. The threshold current 
density amplitude is of the order of the dc current density switching the magnetic 
junction. 
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1. Introduction 

There is unremitting attention to the behaviour of magnetic junctions under spin- 
polarized current flowing. It is not surprising because of a number of interesting 
phenomena which have been observed, such as magnetic conflguration switching [1], 
spin wave generation [2], etc. The effects can occur on a nanosize scale, since their 
characteristic lengths are the exchange and spin diffusion ones with typical values of the 
order of 10 nm. This allows to use such effects for high-density information recording by 
electric current, unattainable for switching magnetic elements by magnetic fleld alone. 

The current-driven switching of magnetic junctions is accompanied often with 
magnetization oscillations and the other high-frequency effects (see, e.g., [2]^[5]). In this 
connection, an interesting problem arises, namely, effect of spin-polarized alternating 
current on magnetic junctions. 

In this work, we consider an effect of alternating (high-frequency) current on the 
magnetic junction conflguration. When the parametric resonance conditions fulflll, both 
collinear conflgurations, parallel and antiparallel ones, can become unstable. It should 
be noted that the parametric resonance in magnetics was studied in many works (see, 
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e.g., PI). However, the parametric resonance was considered there which was excited 
by high-frequency external magnetic field, i.e. the nonlocal Ampere field. In our case, 
the spin-polarized current interaction with magnetic lattice is of exchange nature, so 
that it is localized in the above-mentioned small range. As will be shown below, the 
current-induced parametric resonance has additional features. Thus, the main effect 
takes place at the precession frequency, not the doubled one. The instability is possible, 
too, beyond the parametric resonance conditions. It appears that only the spin-injection 
mechanism [TIIH], not the spin-torque transfer (STT) mechanism P,[in], contributes to 
the effects in consideration. 



2. Model considered and main equations 



We consider a conventional spin-valve model consisted of a pinned ferromagnetic layer 
(layer 1), thin spacer layer, ferromagnetic free layer (layer 2) and nonmagnetic layer 
(layer 3) closing electric circuit. The alternating current flows perpendicular to the 
layer planes (CPP mode). We investigate stability of coUinear (parallel or antiparallel) 
relative orientation of the pinned and free layers against small magnetization fluctuations 
under alternating current flowing with density 

j{t) = jo cos nt. (1) 

The free layer is assumed to be thin compared to spin diffusion length and domain 
wall thickness, so that the macrospin approximation is applicable [11]. In this 
approximation, the fluctuations are described by the modifled Landau-Lifshitz-Gilbert 
(LLG) equation [TT] 
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Here the following notations are used: M = M/\M\ is unit vector along the free layer 
magnetization M, Mi is the same for the pinned layer, H is external magnetic fleld, 
i^a is anisotropy field, n is unit vector along the anisotropy axis, is demagnetization 
field, L is the free layer thickness, a is the magnetization diffusion constant, 7 is 
gyromagnetic ratio, k is the Gilbert damping constant. The parameters p and k describe 
the spin-polarized current effect on the free layer magnetic lattice due to spin-injection 
mechanism [71 E] and STT mechanism [9l [10] , respectively. In coUinear configuration 
with (^Mi ■ M^ = ±1, these parameters take the following forms: 
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where e is electron charge, /ib is the Bohr magneton, Qi is the conductivity spin 
polarization in layer 1, A = L/Z <^ 1, Z being the spin diffusion length in the free layer, 
r is spin relaxation time in the free layer, a is the dimensionless sd exchange interaction 
constant in the free layer, b = (aiMiri)/(aMr) ratio describes the pinned layer 
contribution (see [12] for more details), z/ = Z1/Z2, v* = Xu + (Z1/Z3); Zi {i = 1, 2, 3) 
being the spin resistances [13] 

= adi-Qir 

ai is conductivity of the i-th layer. The upper and lower signs in ^ refer to the parallel 
and antiparallel configurations, respectively. 

The formulae (I2l)-(jll) have been derived on the assumption of direct current 
flowing. However, when the alternating current frequency VL is comparable with 
the magnetization precession frequency, the conduction electrons can follow the 
magnetization oscillations, so that all the conditions fulfill under which the formulae are 
valid. Therefore, we may substitute the alternating current density ([1]) with frequency 
Vt for j in formulae ([3]) and (jlj). Then the parameters p and /c in ([2]) will have time 
dependence of the form p'^'^\t) = p'g^^l cosf2t| and k{t) = kocos^t, respectively. 

Note that the contribution of the spin-injection mechanism is proportional to the 
absolute value of the current density, so that it is the same for forward (1 — 2 — 3) 
and backward (3 ^ 2 ^ 1) currents, in contrast with the contribution of the STT 
mechanism. This leads to different spectra of p^'^^t) and k{t) functions: 

=p[,^)|cosf)t| 

... /2 4 4 \ 
= - H cos 2nt cos4fit + . . . , (6) 



,7r Stt IStt 
i.e., only the even harmonics of Q frequency (including dc component) are presented in 
the p^'^^t) spectrum, while the k{t) function spectrum consists of the single frequency f2. 
Therefore, two mechanisms do not interplay in the lowest-order resonance phenomena 
and they may be considered separately. We begin with the spin- injection mechanism. 

Let X axis be directed along the current, yz plane be parallel to the layer planes, 
the free layer occupy the range < x < L, vectors H, n and Mi have the following 
components: H = {0, 0, i^^}, n = {0, 0, 1}, Mi = {0, 0, 1}. We investigate stability 
of equilibrium collinear configurations = ±1 against the free layer magnetization 
small fluctuations M^, My. The LLG equation components linearized in the fluctuations 
taking the spin-injection mechanism into account only, take the form 

^ + '^^''^ + + 

. 0, (7) 



Alternating current driven instability in magnetic junctions 



ap^^\t) 



Mr = 0. 



(8) 



The periodic time dependence of the coefficients of the last terms in the left-hand side 
of the equations leads to possibility of parametric resonance. 

3. Parametric resonance 



It is well known [Tl], that the parametric resonance occurs when the parameter 
modulation frequency is close to 2ujQ/n, where uq is eigenfrequency of the oscillatory 
system, n = 1, 2, 3, . . .. If the modulation coefficient is small, the parametric instability 
range narrows and the instability threshold rises with increasing the resonance order n. 

In accordance with iQ, we consider parametric excitation at frequency 2Q with the 
ffist two terms taking into account in the right-hand side of ([6]). We assume the damping 
constant k, to be small and neglect it for the time. Taking the Fourier transforms 
of (CI), (IHl) with respect to time, we have 
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is a quantity proportional to the current density amplitude with dimension of frequency. 
Usually, a condition is fulfilled 

47rM> H^, 3eW/7, (12) 

that is assumed below. 

If we make the substitution u ^ u ± 2Q in (ITOl) . (fTTj) . the equations are obtained 
where terms with Mx^y{Lj±2Q) stand in the left-hand sides and the terms with Mx^y^iv) 
and Mx^y{uj ± 4^2) with e*-^-* coefficient in the right-hand sides. 

If frequency e^^-* (the coupling parameter) is small compared to the other 
characteristic frequencies {uq, Q), the coupling with M^^y^uj ± 4f2) can be neglected, 
because it leads to higher-order corrections in e^^^. As a result, we obtain a closed 
system of equations for Mx,y{^^) and Mx,y{uJ ± 20). Equating the determinant of the 
system to zero, we get the dispersion equation 
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(13) 
(14) 



is the eigenfrequency neglecting the terms quadratic in e'-^-', 

A{uj) =ujI-uj^. (15) 

Note, that the dc component ofp^'^^t) (see ([6])) renormahzes the system eigenfrequency. 

Near the parametric resonance, we have u ^ ujq ^ Q, so that we may retain only the 
summands with resonant denominator A{uj — 2il) in the terms with e*^^-* and replace u 
and Q with ujq everywhere but that denominator. This leads to the parametric resonance 
equation 

A(cj)A(cj - 2^) = {AiT-fMe^^'^f. (16) 
Let u = ujq + u, Q = ujq + 6, where |z/|, \6\ ^ ujq. The equation for u takes the form 
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It is seen from ( IT8|) that the parametric instability takes place at 
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In presence of dissipation (k 7^ 0), damping takes place with decrement (see, e.g., 0) 
|(Imcj)d| = 27rK7M. (21) 

If K <^ 1, the parametric instability threshold may be estimated from the condition 

(Imcj)p > |(Imcj)d| (22) 

at zero resonance detuning {6 = 0). 

Equations fl20l) - fl22|) give the following condition for the parametric instability: 

e^^) > kcoq. (23) 

The right-hand side of this inequality contains e*^^\ too. Therefore, the inequality is 
to be resolved with respect to e*-^-*. As a result, the following instability threshold is 
obtained for the collinear configurations 
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In contrast with the case of direct current, both colhnear configurations can become 
unstable, but the corresponding thresholds are different. At e*^") > e[^\ e^^^ < e[^^ the 
switching from unstable antiparallel configuration to stable parallel one is possible, while 
at e^"*"-* > e[^^ both configurations are unstable. At i^a =t -^2 ^ Ottk^M the instability 
threshold for both configurations takes the form 

eiJ) = 2K^[7rM{H., ± H,)]'/^ = kJ^\ (25) 

where uJq^^ = 2'y[TTM(H^±Hz)Y^'^ is the eigenfrequency in absence of the electric current. 
If H^- <t: gvTK^M, we have 

4"^ = lliH, - H.) (26) 

for the antiparallel configuration. It is seen from (12^ . that the instability threshold can 
be lowered considerably with external magnetic field close to, but slightly lower than the 
anisotropy field. Note, that participation of the magnetic field does not prevent locality 
of the effect, because the magnetic field cannot do switching alone, without the current. 

Let us compare the alternating current density amplitude jth corresponding to the 
parametric instability threshold with the the direct current density jdc leading to the 
switching antiparallel orientation to parallel one in absence of the external magnetic 
field. The dc threshold corresponds to the condition [12] ap^~^/L > '~^H.^, where p*^") is 
determined by ([3]). In the parametric resonance case, a'p\^ jL = ?)t:k'-){t:MH^^/'^ ^ so 
that 

,,^J^f. (27) 



jdc \Hg,J 

At typical values of the parameters {M/ H^, ~ 10, k = 3 x 10^^), this ratio is of the order 
of 1. At lower damping constant, the parametric instability threshold will be smaller 
than the dc threshold. 



4. Non-resonance instability 

The instability of the antiparallel configuration under alternating current fiowing is 
possible also when the parametric resonance condition Vt ^ ujq does not fulfill. It follows 
from (|T4l) that the eigenfrequency ojq becomes imaginary at e*^^' > 'j{Hi, — Hz) /3 because 
of the contribution of dc component in the spectrum of p^'^^t) function, i.e., such a 
component of the sd exchange field. The cancelling of the eigenfrequency corresponds 
to an orientational phase transition similar to that under dc injection current [I2]. The 
threshold amplitude of the alternating current is 7r/2 times as much as the corresponding 
dc threshold. Note, that the threshold is higher then the parametric instability threshold 
e[^\ so that the parametric instability develops first under fulfilled parametric resonance 
conditions. 
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5. Is similcir effect possible due to the spin-torque transfer mechanism? 

Since the k{t) function describing the STT mechanism contribution has a single-mode 
spectrum, the lowest order of the parametric resonance corresponds to Q ~ 2u!o 
condition. However, if all the previous calculations are carried out for p^'^\t) — 
0, k{t) 7^ 0, the contributions from k{t) mutually cancel in resonance approximation 
{Q, 2u; 2u;o). Therefore, a negative answer should be given to the question stated 
above. This is due to the fact, that the STT, in contrast with the spin injection, 
modifies the damping, not the eigenfrequency. As is known, the parametric resonance 
is not possible in such a situation. 

Such a difference between two mechanisms can be used for separating their 
contributions in various cases. 

6. Conclusion 

The analysis carried out shows that the coUinear configurations of a magnetic junction 
can become unstable under parametric resonance conditions. The instability threshold 
is of the same order of magnitude as the corresponding threshold under direct current 
flowing through the magnetic junction. Depending on the alternating current density 
amplitude, the antiparallel configuration can be switched to parallel one or both coUinear 
configurations can be unstable with growing magnetic fluctuations. To elucidate the 
resulting state which the system considered comes to, further investigations, both 
theoretical and experimental, are needed. 
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